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1 Introduction 



The study of eigenfunctions of quantum systems, in particular their dependence on the classical 
dynamics, has attracted a lot of attention. A prominent class of examples is provided by two- 
dimensional billiard systems, which are classically given by the free motion of a particle inside 
some domain with elastic reflections at the boundary. The corresponding quantum system is 
described by the Helmholtz equation inside a compact domain f2 C (in units h = 1 = 2m), 

AV^„(a;) + kl^nix) = , xeVl , (1) 

with (for example) Dirichlet boundary conditions 

ipnix) = , xedQ , (2) 

where the eigenfunctions ijjn{x) are in L^(f2). Assuming that the eigenvalues k"^ are ordered 
with increasing value, the semiclassical limit corresponds to n —>■ oo. A detailed knowledge of 
the behaviour of the eigenvalues fc^ and the structure of eigenstates is relevant for applications, 
for example microwave cavities or mesoscopic systems (see e.g. |1| and references therein). 

For the description of the phase space structure of quantum systems usually the Wigner 
function [2| or Husimi function 3J are used. However, for a system with d degrees of freedom 
these are 2(i-dimensional functions, which are difficult to visualize for d > 1. Therefore, one 
usually considers the position representation, or the momentum representation jl], or sections 
through the Wigner or Husimi function, see e.g. 

Another approach is the use of representations on the billiard boundary, acting as a global 
Poincare section. In the literature one can find several variants for these representations, see 
e.g. inilZlE]- The reason is, as emphasized in [7], that there is no natural definition of a scalar 
product for functions on the billiard boundary. This raises the question whether one of these 
definitions has advantages over the others, which will be addressed in the following. 

The representation of eigenstates on the Poincare section plays an important role in several 
applications. For example, it is used to define scar measures [HI El, or to study conductance 
fiuctuations, see jTU] and references therein. Furthermore, these representations are used to 
determine the coupling of leads in open systems JT]. Another important application is the 
detection of regions where eigenstates localize, see e.g. [T21IIS1III] (for an alternative approach 
based on the scattering approach see OCS]). Representations of eigenstates on the Poincare 
section have also been useful to understand the behaviour of optical microresonators, see e.g. [in] 
and references therein. More generally, the approach is not just applicable for billiard systems 
but it is also useful for Poincare sections arising from Bogomolny's transfer operator approach 

El- 

In this paper we first compare two different definitions for the Poincare Husimi representa- 
tion, discuss their properties (section |2)) and based on this we select one particular definition 
for the following. In section [31 we derive the behaviour of these Poincare Husimi functions when 
averaged over several energies. In section [^ we establish a relation between the well-known 
Husimi function in phase space and the Poincare Husimi function on the billiard boundary. 
This allows for a direct physical interpretation of the Poincare Husimi functions. Moreover, for 
ergodic systems a quantum ergodicity theorem for the Poincare Husimi functions is shown. 
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2 Husimi representation on the boundary 



Let us first recall the definition and some properties of Husimi functions in phase space. For a 
solution ipn of the Helmholtz equation ((H) with energy E = the Husimi function H^{p, q) 
is given by its projection onto a coherent state, i.e. 

^f(P,q):= i^)' Mp,.)U^^M' ■ (3) 

Here 

{^i,^2)n:= I M<l)Mq)d'q (4) 



is the scalar product in Q, and ipi denotes the complex conjugate of ipi. 
The coherent states are defined as 

V^(p,q),fc(=^) := (^)'^'(detImi?)i/2g^'=[(p,.-q)+i(.-q,B(.-9))] ^ 

where (p, g) G x denotes the point in phase space around which the coherent state 
is localized, and 5 is a symmetric complex 2x2 matrix which determines the shape of the 
coherent state. For the conventional coherent states one has i? = i ( q 5 ) an in general one 
has the condition Imi? > 0, i.e. {v,lmBv) > for all v G -R^\{0}. Notice, that because 
the variance of the coherent states is proportional to k, all Husimi functions are concentrated 
around the energy shell |pp = 1 (and not around |pp = k'^). By this it is possible to compare 
Husimi functions with different energies k"^, and for example consider their mean, see ((Tj) below. 

Such Husimi functions can be interpreted as probability distributions on phase space, be- 
cause they satisfy the relation 

Ai;^)n = j j a{p, q)H^{p, q) d^pd^g + 0{k-^) (6) 

R2 R2 

where a(p, q) is a function on phase space and A its quantization. Moreover, the average of all 
Husimi functions H^{p, q) up to some energy k"^ = E converges for A; ^ 00 to the normalized 
invariant measure on the energy shell. 

Here N{k) denotes the spectral staircase function, A^(A;) := < k}, and Xn{Q) is the 

characteristic function on Q. The mean behaviour (|7j) is similar to the mean behaviour of 
the spectral staircase function, which is given by the Weyl formula, i.e. for k ^ 00 one has 
N{k) ~ ^k'^ where A is the area of the billiard. A similar asymptotic behaviour can be derived 
for the mean of normal derivative functions, see JH] for a detailed discussion. 

For billiards an extremely useful approach for describing the dynamics is the use of a 
Poincare section V together with the corresponding Poincare mapping P. Usually the sec- 
tion V := {{q,p) I q G [0, vol(9fi)], p G [—1, 1]} is parameterized by the arc-length coordinate 
q along the boundary dQ of the billiard and the projection p of the (unit) momentum p af- 
ter the reflection on the tangent t{q), i.e. p = {p,t{q)). By this the billiard flow induces an 
area-preserving map P : V ^ V where the invariant measure is given by d/x = dg dp. 
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In order to have the advantages of such a reduced representation in quantum mechanics as 
well, one is interested in a Husimi representation hn{q,p) on the Poincare section V which is 
associated with an eigenstate ipn- Such a Poincare Husimi function should have similar pro- 
perties as the ones expressed by relations © and ((Tj) for the Husimi functions in phase space, 
and our aim is to study to what extent this is possible. More precisely, one would like that for 
the Husimi function on the billiard boundary a spectral average 

T<k{q,p) ■■= j;^^hn{q,p) (8) 

tends to the invariant measure on P as — > cx), in the same way as in ((Zj). 

The Husimi representation on the billiard boundary is usually defined using the normal 
derivative of the eigenfunction (hereafter called the boundary function) 

M„(s) := (n(s),VV^„(a;(s))) , (9) 

where x{s) is a point on the boundary (9fi, parameterized by the arc-length s and n(s) de- 
notes the outer normal unit vector to dVt at x{s). The boundary functions are a natural 
starting point for defining a Husimi representation because they determine the eigenfunctions 
uniquely, see (j^ . Thus the boundary functions form a reduced representation of the system. 
If an eigenfunction ipn is normalized, then the corresponding boundary function Un fulfils the 
normalization condition 



\u^{s)\''{n{s),x{s))ds = kl . (10) 



do. 



For alternative derivations of (fTUj) and more general boundary conditions see [221 123 • Notice 
that while the integrand depends on the chosen origin for the vector £c(s), the integral is 
independent of this choice. 

Starting from the boundary function a Husimi function on the Poincare section can be 
defined by a projection onto a coherent state. There are different possibilities to define coherent 
states on the boundary of a billiard. A natural choice is the periodization of the usual one- 
dimensional coherent states. 



where (g,p) G dVt x M, and L denotes the length of the boundary. The parameter 6 G C, 
Im 6 > 0, determines the shape of the coherent state. Then for an eigenstate ipn with boundary 
function Un a Husimi function on the Poincare section V (or more precisely, on the cylindric 
phase space dVL x M) can be defined as [7j 



hn{q,p) 



2 

1 



an 



2nk. 

The completeness relation for the coherent states gives 



(12) 



J J hn{q,p)dpdq = ^ j |u„(s)pds, (13) 



an R an 
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Figure 1: Examples of eigenstates ipnio)-, shown to the left, and to the right their 
Poincare Husimi functions hn{q,p)- In a) an eigenstate (n = 1952) localizing around 
a regular orbit for the limagon billiard at e = 0.3 is shown. In b) and c) two 
eigenstates for the cardioid billiard are shown (n = 1817 and n = 1277). 



so in view of relation ()10|) the Poincare Husimi function hn{q,p) will in general not be nor- 
malized. This can be fixed by dividing hn{q,p) by the factor p- J|M„(s)pds, as was done 



for instance in [121 II2]- But later on we will see that it is more natural to work with the 
non-normalized Husimi functions p2|) . 

A different Poincare representation has been proposed in |H], 



hn{q,p) 



I C(g,p),fc„(s) Unis) {n{s),x{s)) 
an 



(14) 



an 
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P -1.5° 



b) 




Figure 2: Plot of 'Hk{q,p) for k = 125 using the first 2000 eigenstates in the hmagon 
bilhard of odd symmetry at e = 0.3. In a) the result for Hkiqyp) using definition 
(fT^ for hn{q,p) is shown and in b) a corresponding Hk{q,p) using definition (fT^ is 
displayed. In addition to the symmetry related dips at {q,p) = (0,0) and (L/2,0) 
one clearly sees the variation in p-direction in both cases and in b) we moreover 
observe a variation in q. 
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where the inclusion of the factor {h{s) , x{s)) is motivated by its appearance in the nor- 
mahzation condition In order to compare the two definitions, we use the fact that 

for large k the coherent state becomes more and more concentrated around s = q and so 
(n(s),a;(s))cj^p)^^js) ~ (n(g), a3(g)) cj^p^^^js). This leads to the relation 

hn{q,p) ^ {n{q),x{q)) hn{q,p) (15) 

between the two definitions for Husimi functions. 

Let us first illustrate the behaviour of the Husimi representation given by (jl2|l . As a concrete 
example we consider a member of the family of limagon billiards introduced by Robnik [24, 25], 
whose boundary is given in polar coordinates by p{ip) = 1 + e cos{ip) where e G [0, 1] is the 
family parameter. At e = 0.3 the billiard has a mixed phase space (see figure 1 in J2]) and at 
£ = 1 it turns into the fully chaotic (i.e. ergodic, mixing, . . . ) cardioid billiard. Because of the 
symmetry of the billiard we consider the half limagon billiard with Dirichlet boundary conditions 
everywhere. Fig. ^ shows a comparison of eigenstates ipn{<l) with their Husimi representations 
hn{q,p) as grey-scale plots with black corresponding to large values. For the computations 
b := ia~^ = i was chosen. In a) an eigenstate which is localized around a stable periodic orbit 
with period three is shown which is clearly reflected in its Poincare Husimi function to the 
right. The symmetry hn{q,p) = hn{q, —p) is due to the time-reversal symmetry of the system 
and the symmetry hn{q,p) = hn{L — q,p) stems from the reflection symmetry of the system. 
The plots in fig. ^) and c) are at e = 1.0, i.e. for the cardioid billiard. The eigenstate shown 
in b) is localized around an unstable periodic orbit of period two which is also nicely seen in 
the prominent peaks for the corresponding Poincare Husimi function. In c) an irregular state 
in the cardioid billiard is displayed which is spread out over the full billiard and also hn{q,p) 
does not show any prominent localization. 

Now we turn to a comparison of the two Poincare Husimi representations given by p2j] 
and (fT^ . In figure |2 a plot of 'Hk{qiP) is shown where k = 125.27. . . is chosen such that the 
first 2000 states are taken into account. Both definitions, equations (fT^ and (HD), lead to a 
similar non-uniform behaviour of 'Hk{qiP) in P direction. We will discuss this behaviour in 
more detail in the following section. In addition we observe that Ti.k{q,p) has a minimum at 
{q,p) = (0,0) and (g,p) = (±£/2,0) which is due to the desymmetrization. Figure|2b) shows 
a plot of 'Hk{q,p) which is defined as 'Hk{.'i,p)i but instead of hn{q,p) the functions hn{q,p) 
are used, see definition ()14|). In this case we observe in addition a clear variation in q. The 
reason for this is the factor {n{q), x{q)) as explained by relation ()15|) . Another important point 
is that the definition ()14|1 depends on the chosen origin as the factor {fi{q),x{q)) does, and 
therefore the integrals in equation (fT^ are not invariant under a shift of the origin. Because 
of the variation of hn{q,p) in q and the dependence on the origin we prefer the definition (fT2)) 
and will use this exclusively in the following. 



3 Mean behaviour of boundary Husimi functions 

In this section we determine the asymptotic behaviour of the mean 'Hk{qiP) of the boundary 
Husimi functions for large energies. To this end we will use the methods from our previous 
work [in]. Let us introduce 

f{k, s, s') := J2 ''"^'|!"^'^^ p(fc - fcn) , (16) 



7 



where p is a smooth function whose Fourier transform is supported in a neighbourhood [—rj, rj] 
with ry smaller than the length of the shortest periodic orbit of the billiard flow. The function 
g^{k, s, s') was studied in [12] and an asymptotic expansion was derived. Its leading term reads 



27r 



27r2 



(n(s),e(^))(n(s'),e((^))e''<"(^)-"(^')'^('^» dip (1 + 0(A;-^)) , (17) 



where x{s) denotes the position vector on the boundary at point s, n(s) denotes the outer unit 
normal vector to the boundary at s and e{ip) = (cos 93, sin is the unit vector in direction cp. 
Multiplying (fTTjl with c^^p-j ^(s) and c^^p-j ^(s') and integrating over s and s' leads to 



^p{k - kn)hn{q,p) 



neN 



k^ 
47^3 



2n 



(n(s),e(^))e''=<-W'^(^»4,),,(s)ds 



d^ (l + 0(A;-i)) . 



(18) 



an 

The s-integral can be computed by the method of stationary phase, 



(n(.),e(^))e^'^<-W'^(-))4^),,(.)d. 



an 



- j (Im 6)^/^ j (n(s), e{ip)) eiM(-(«),eM>-P(«-<?)-5{--9)'] ds 
—00 

^y/'i^jl^ (n(g), e(v.)) e^MM.),e(^)>+i(P-(t(.),e(^)»^] + Qik-'/')) 



(19) 



with 



6 = 6 + K(g)(n(g),e(v9)) , 
where K{q) is the curvature of the boundary at q. Inserting this result we obtain 

^p{k - kn)hn{q,p) 



(20) 



2n 



2P Mtt\i f (Imb) 



(27r)3 \ k J 



n/2 



1^1 



(n(g), e(^)) r e-'^tF^^-^*^^)'^^^)^)^ d^ (1 + 0(A:-V2)) 



(21) 



and for \p\ < 1 the (y9-integral can again be solved by the method of stationary phase (notice 
that there are two stationary points) which yields 



k 



Tpik - kn)K{q.p) = -^1^ (1 + Oik-'/')) . 



(22) 



neN 



By integrating this equation, and using a delta sequence for p as in proofs of the Weyl formula 
(see e.g., |2E]), we finally obtain 



T-(-k{q,p) = 5Z ^'^(^'P) = ^ Vl-P^ + C'(r 



1/2N 



(23) 
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In the derivation of ()22|) from ()21|) we have assumed that \p\ < 1 because then the stationary 
points are non-degenerate. For \p\ > 1 the stationary points become complex and the integral 
is exponentially decreasing for k oo. 

Previously, such a -y/I — behaviour appeared in the context of Fredholm methods for 
Poincare Husimi functions |2Z1 and was also obtained in connection with the inverse participa- 
tion ratio p. 

Next we want to derive a uniform approximation which describes the mean behaviour of 
the Husimi functions near \p\ = 1 and the crossover from the regime \p\ < 1 to the exponential 
decrease for \p\ > 1. We will study the case p ~ 1, the case p — 1 is completely analogous. 
Let (fo be the angle corresponding to the direction of t{q) and expanding the amplitude and 
phase function in (^1]) around ipo leads to 

oo 

, \u / \ 4P /47r\l/2 f (Im6)l/2 2)2 r^n -l/2^^ 

2_^p{k-K)K{q,p) = j—^i^—) / ' ~ f^e ^ > d^{l + 0{k '^^)) 

oo 

/T 7 \ 1 /9 ^ 



(24) 



{27r)^\kJ J \b\ 



(\ 1/4 oo 
K\ /../v^<'-"-*d.(i + o(r./^)) 
Imo / J 

where D_3/2{x) denotes the parabolic cylinder function and we have used one of the standard 
integral representations, see e.g. j2Hl- 

This result was derived under the assumption p ^ 1 such that (p^ — 1) ^ 2{p — 1). Sub- 
stituting [p — 1) by (p^ — l)/2 allows to combine the results for the different p-regions in one 
formula 

J2pik - k^)K{q,p) = A F,{p) (1 + 0{k-'/')) , (25) 

neN 

where 

/ ~ \ 1/4 



FJp) = ———eW^'^-' \^\ ^-3/2 ^ (p' - 1) • (26) 

2(2A;)V4 llmfe/ V 2V2|5| ^ ^ ^ 



For IpI < 1 one has Fk{p) = a/1 — p^ + 0(A; ^), since Z)_3/2(x) ~ 2^/^|x|"^/^e^ for x — oo. 
Recall, that b is defined in equation ()20|). In figure El we compare the expression with 
Im6 = 1 for different values of k. It is clearly visible that the asymptotic result is reached 
slowly with increasing k. 

Integrating (j^ . analogous to the transition from to ^I^ . one can compare the uni- 
formized mean behaviour with the numerical result. In figured a section of 'Hk{q,p) at q = 3.0 
is shown for k = 125, compare with figure |2ti') • The remaining differences are due to higher 
order corrections. 

In the derivation of the results ()22|1 and ()25j) we have implicitly assumed that the boundary 
of Q is sufficiently smooth, because only then we can use the stationary phase formula. But it 
is easy to extend the results to the case that the boundary is only piecewise smooth. Since we 
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0.0 0.2 0.4 0.6 0.8 1.0 1.2 p 1.4 



Figure 3: Comparison of the uniformized asymptotic behaviour Fk{p), see fl2(jj) . with 

= 1 and for k = 10, 30, 500. The asymptotic semi-circle behaviour is reached 
slowly. 

multiply in (|18|) by a coherent state centered in q, all the following computations remain valid 
if q is in the smooth part of the boundary, since the contributions from the singular points are 
exponentially suppressed then. So it could only happen that some additional mass sits at the 




Figure 4: The full curve shows a section of Ti.k{q,p) at g = 3.0 with k = 125 for the 
desymmetrized limagon billiard, see figure |2t), and the second line is the uniformized 
mean behaviour. The remaining deviations are caused by higher order corrections. 
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singular points of the boundary, i.e. that we have 



hm -T—TT-r 

fc^oo N{k) 



Yl hn{q,p) = ^ + fisip, q) , (27) 



where fis{p, q) dpdg is a measure supported on the singular part of the boundary. Since hn{q,p) 
is positive, and \/\ — p^ is absolutely continuous with respect to the Lebesgue measure, we 
have /i5 > 0. ^ Now the completeness relation gives limfc„_^oo | // {"f^iq), x{q))hn(yq,p) dpdq = 1 
and i JJ {niq), x{q))-^ ^Jl — dpdg = 1 and therefore 

1 

j J {niq),x{q))fis{p,q)(iqdp = o . (28) 
-1 an 

But for a star-shaped billiard one can choose the origin of the coordinate system such that 
{fi{q),x{q)) > for all q G dQ, and so Hs = 0. Therefore ()22j) and ()25|1 remain true for 
star-shaped billiards with piecewise smooth boundary with the only possible modification that 
the error term might decay more slowly at the singular points of the boundary. 

4 From Husimi functions in phase space to Husimi func- 
tions on the boundary 

In this section we derive a direct relation between the Husimi function in phase space and the 
one on the Poincare section, as given by eq. p2|l . By this we obtain a physical interpretation 
of the Poincare Husimi representation. For the calculations in this section we have to assume 
that the billiard domain Q is convex. Let he a solution of the Helmholtz equation (^J) in 
which satisfies Dirichlet boundary condition on dfl. Any such function can be represented as 

^Ij{x) = - j Gk{x-x{s))u{s) ds (29) 
an 

where Gk{x — y) is a free Green function and u{s) is the normal derivative of ip on the boundary. 

Let ipz be a coherent state © centered at z = (p, q) G T*M^, for reasons of simplicity we 
restrict ourselves to the case of a non-squeezed sjTiimetrical state, i.e. B = i (q 5), and omit 
the index B in the following. We want to compute the overlap {ip, ipz)-, given by 

(V', i'.)n = - j {Gk{- - x{s)),i,z)nu{s) ds , (30) 
an 

and we observe that 

{Gu{--x{s)),^ljz)n = Gl^.{x{s)) (31) 

where 

Gfc = lim— — ^— (32) 

^To see this, assume that has a negative part, then there exist a point zq = {pq, qq) and constants Eq, C > 
such that //|^_^g|<g/isd2: < -C for aU s < So- This imphes hme_+o //|^_^„|<g linife^oo iv(fe) Z]fc„<fc ^n(^) < 
— C, which contradicts the positivity of the /i„(z). 
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is the resolvent operator, whose kernel is the Green function. From equation pip we see that 
the function G^^ipz is restricted to the billiard boundary. For the resolvent operator we use the 
integral representation 







Gt = i y e'^* U(t) dt (33) 

— oo 

where U(t) = e^*^ is the free time evolution operator with 1/k playing the role of h, and 
inserting equation into we obtain 



= ^ J e"^'Vit)i;,ixis)) dt . (34) 

— oo 

But the free time evolution of a coherent state centered in z is well known (see e.g. [^ 1^0]) 
to give again a coherent state, centered around the image of z under the classical flow and with 
transformed variance, 

V{t)^z{x) = e'^P'* (^^^ giM(P.-'?{t))+5(Ti2u){--'?{*))'] , (35) 

with q{t) = q + 2tp. Therefore, G\ipz{x) has the structure of a Gaussian beam emanating 
from the point q in direction p backwards in time. If we introduce a new coordinate system 
X = (a;||,a;_L) centered at q with x\\ parallel to p and x± perpendicular to p, we obtain by a 
stationary phase approximation that for x± and 1 — |p| small (i.e. near the energy shell) 

^i"'-*-* = r.,nl ^,„ e""-"^'^-'^''-"-"'' (1 + 0(*-'/^)) (36) 

V2A;(1 + 

holds, where we have assumed that x\\ < 0. For a;|| ~ and X|| > the integral leads to an error 
function which describes the transition from the exponentially decaying regime with x\\ > to 
the regime X|| < in ()Hfi|l . For |p| = 1 the result reads 



ik[x 



^2 1 



n^2(i+ix||)''±j 



|p|=i V2k(l + ixuV/^ 

^ " (37) 



where erfc(2;) denotes the complementary error function, and the absolute value of this expres- 
sion is shown in figure El 

Next we want to evaluate this expression on the boundary. To this end, let x{q) be the 
point of intersection between the boundary and the line from q in direction —p. (Here we need 
the assumption that the billiard domain Q is convex, in order that there is only one such point.) 
Then we obtain with x{s) = x{q) + i{q){s — q) — 'h{q){s — qY + 0{{s — q)^) that 



x^ = {l-pY\s-q) + 0{{s-q)') (39) 



X|| = \q- xiq)\ +p{s~q)-^{l- p'f'^s - qf + 0((s - qf) (38) 
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Figure 5: Illustration of a Gaussian beam as given by pTj) inside the limagon billiard 
at e = 0.3. 



where p := (p, i) G [—1, 1]. Inserting these expressions in (jHBjl gives 



■ 1/4 1 



where c^gp^^(s) is a coherent state on the boundary, as defined in (fTT|) . with variance b = 

- - P^)^^^' = (|g-"xfa)|^+i)i/4 - Notice that although we started with a 

symmetric coherent state in the interior, the projected coherent state on the boundary is no 
longer symmetric and has a non-trivial squeezing parameter b which depends on the position 
of the original state, the angle of intersection of the ray in direction —p with the boundary and 
the curvature of the boundary. 

If we insert the expression ()40|) into (jHUj) we obtain a semiclassical relation between the 
projection of an eigenstate onto a coherent state in the interior and the projection of the 
normal derivative on the boundary onto a coherent state on the boundary. 



i7ri/4 1 



gifc„|q-a,(g)l+ieg-i^(l-|p|)2 



(41) 
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In turn from this we obtain the central result of this section, a direct relation between the 
corresponding Husimi functions 



Hnip, q) = -\P\)\ (1 + OiK'/')) , (42) 

with 

1/2 

e-fcnd-IPl)^ . (43) 



Let us first discuss the meaning of the individual terms on the right hand side of equation 
(B^ . The function Sk„{l — \p\) is a delta sequence for /c„ — >• oo, and describes the localization 
of Hn{p, q) around the energy shell. The factor 1/ a/1 — comes from the projection of the 
Gaussian beam to the plane tangent to the boundary, see figure El 

As in the previous section we have assumed that the boundary is smooth. But by the 
localization of the coherent states the results can be again extended to the case that the 
boundary is piecewise smooth, then remains valid if q is not a singular point of the 
boundary. 

The direct connection between the Husimi function in the interior and the one on the 
boundary, given by equation (jl^ . allows to derive interesting relations between the two Husimi 
functions and can be used to give a direct physical interpretation of the Husimi function on the 
boundary. From equation © together with relation (fl^ we obtain 



1 

(V^„,AV^„)n= I [ -^^^)M=(a)(g,p)/(g,p)dgdp + 0(fc-V2) , (44) 
4-\/l — 



-1 an 



where l{q,p) denotes the length of a ray emanating from q{q) G dQ in the direction determined 
by p until it hits the boundary again. Furthermore, 

Kq,p) 

{a){q,p) := -j^^—^ j a{q{q) + te{q,p), e{q,p)) dt , (45) 


is the mean value of the classical observable between two bounces, where e{q,p) denotes the 
unit vector at q{q) in direction p. A relation of the same type as has been obtained recently 
by different methods in [3T] for certain localized functions on the boundary. 
We conclude from relation ()44|) that 

1 K{q,p) , . 

Wg,p):=4^7= (46) 

is a reduction of the probability density defined by the Husimi function on the whole phase 
space to the boundary. So if one wants a proper representation of eigenf unctions on the Poincare 
section which is an approximate probability density, and whose general properties are indepen- 
dent of the billiard shape, then ()4fi|l seems to be the best choice. Of course a drawback of the 
function ()4fj|l is the singularity of 1/ -^/l — at p = ±1 which is relevant at any finite energy. 
So for numerical computations the definition (fT^ is more suitable and the importance of ()46j] 
lies in the physical interpretation. 
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In particular, relation ()44|) implies an asymptotic normalization condition on i)n{Q,p), 

1 

Uq,P) Kq,P) dgdp = 1 + Oik-'/') . (47) 
-1 an 

Since l{q,p)dqdp is the phase space volume in the energy shell corresponding to the volume 
element dqdp of the Poincare section, the factor l{q,p) can be viewed as a normalization which 
makes f)n(g,p) independent of the biUiard shape, i.e. for any V C dfl x [—1,1], we get that 
/x( ^n{<l, p)Kl^ P)^1 is the probability for the particle in the state ipn to be found in the region 
D := U-'V on the energy shell, where the map 11 describes the projection of the domain D to 
the boundary. 

We would like to close this section with some remarks on the implications of quantum 
ergodicity to the behaviour of the Poincare Husimi functions. If the classical billiard flow 
in Q is ergodic, then the quantum ergodicity theorem [321 ES] (see [20] for an introduction) 
tells us that almost all Husimi functions Hn{p,q) tend weakly to 2nvo\{n) • result (jl^ 

then immediately implies that in the semiclassical limit almost all Poincare Husimi functions 
hn{q, p) tend to ^^^^(q,) \/^ ~ P^ the weak sense. So this proves a quantum ergodicity theorem 
for the boundary Husimi functions. Recently related results have been obtained establishing 
quantum ergodicity for observables on the Poincare section Notice that the \/\ — p^ 

behaviour is also visible in the plot of hn{q,p) for the irregular state shown in fig. ^) for the 
ergodic cardioid billiard. 

5 Summary 

Poincare representations of eigenstates play an important role in several areas. However, a 
priori there is no unique way for their definition. In this paper we single out the definition 
given by (fT^ and show that the asymptotic mean behaviour of these Husimi functions is 
proportional to a/1 — p"^. For this asymptotic semi-circle behaviour we in addition derive a 
uniform asymptotic formula. Furthermore we establish a direct relation between the Husimi 
function in phase space and the Poincare Husimi function (jl2j) on the billiard boundary. By 
this a physically meaningful interpretation, see equation (jl^ . of the previously ad-hoc chosen 
definition for the Poincare Husimi function is obtained. Namely, the Poincare Husimi function 
P)„(g,p) can be viewed as a probability density on the Poincare section. For ergodic systems 
our result implies a quantum ergodicity theorem for the Poincare Husimi functions, i.e. almost 
all Poincare Husimi functions become equidistributed with respect to the appropriate measure. 
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